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TECHNICAL NOTE 3653 


SOME EFFECTS OF BLUNTNESS ON BOUNDARY-LAYER TRANSITION 
AND HEAT TRANSFER AT SUPERSONIC SPEEDS 
By W. E. Moeckel 


SUMMARY 

Large downstream movements of transition observed when the leading 
edge of a hollow cylinder or a flat plate is slightly blunted are ex- 
plained in terms of the reduction in” Reynolds number at the outer edge 
of the boundary layer due to the detached shock wave. The magnitude of 
this reduction is computed for cones and wedges for Mach numbers to 20. 
Concurrent changes in "outer-edge" Mach number and temperature are 
found to be in the direction that would increase the stability of the 
laminar boundary layer. 

The hypothesis is made that transition Reynolds number is substan- 
tially unchanged when a sharp leading edge or tip is blunted. This 
hypothesis leads to the conclusion that the downstream movement of tran- 
sition is inversely proportional to the ratio of surface Reynolds munber 
with blunted tip or leading edge to surface Reynolds number with sharp 
tip or leading edge. This conclusion is in good agreement with the 
hollow-cylinder result at Mach 3.1. 

Application of this hypothesis to other Mach numbers yields the 
result that blunting the tip of slender cones or the leading edge of 
thin wedges should produce downstream movements of transition by factors 
ranging from 2 at Mach 3.0 to 30 at Mach 15. The significance of this 
result is discussed with regard to the possible reduction in over-all 
heat-transfer rate and friction drag for aircraft flying at high super- 
sonic speeds. 

Mach number profiles near the surfaces of blunted cones and wedges 
are computed for an assumed shape of the detached shock wave at flight 
Mach numbers to 20. The dissipation and stability of these profiles are 
discussed, and a method is described for estimating the amount of blunt- 
ing required to produce the maximum possible downstream movement of 
transition . 
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INTRODUCTIOIT 

In an investigation of the boundary layer on a hollow cylinder 
alined with the stream direction, Brinich and Diaconis discovered that 
the transition point moved downstream when the leading edge was slightly 
blunted (ref. l) . Similar results were obtained with a flat-plate wing 
in reference 2. A more extensive Investigation of the effects of 
leading-edge geometry on transition (ref. 3) confirmed previous results 
and led directly to the explanation contained herein of the effect of 
blunting on transition. 

When a cone or wedge is blunted slightly (sketch l) , the flow is 



changed in several ways, each of which could have a noticeable effect on 
the transition location. A favorable static -pressure gradient is estab- 
lished near the vertex which could tend to stabilize the laminar layer. 
Downstream of the shoulder, however, the static-pressure gradient is 
adverse (for moderate supersonic speeds) because of the overexpansion 
around the shoulder and subsequent reconrpression to the value corre- 
sponding to the unblunted cone or wedge. The effect of static-pressure 
gradient on transition is therefore Inconclusive. 
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In addition to the static -pressure gradient along the surface, the 
blunting produces a stagnation-pressure gradient normal to the surface. 
This gradient results from the variation in stagnation-pressure loss as 
the detached shock decays from the normal-shock strength at the vertex 
to the strength corresponding to the unhlunted body at some distance 
from the vertex. For in vise id flow, the stagnation pressure along each 
stream-line remains constant downstream of the shock; hence, this gradi- 
ent normal to the surface would persist along the entire length of the 
body. In the constant-static -pressure region a few bluntness diameters 
doimstream of the vertex, the stagnation -pressure gradient results in a 
shear layer whose thickness depends on the size of the blunted portion 
of the body. 

The fact that the entropy gradient produced by strongly curbed 
shock waves might have appreciable effect on the development of the 
boundary layer is' pointed out in references 4 and 5 . Previously, the 
author of the present report had evaluated the shear profiles produced 
by detached shock waves near the surface of blunted flat plates. An 
explanation of the observed movement of treinsltion in terms of these 
shear profiles was therefore sought. 

The interaction of the boundary layer with the shear profile pro- 
duced by a detached shock wave is fundamentally a very difficult analyt- 
ical problem; however, the condition of most interest is one for which 
the interaction of the two profiles is not important. Thus, if the 
shear profile produced by blunting is much thicker than the boundary 
layer, the rate of shear of the former is negligible, compared with that 
of the latter. The boundary layer then develops in a region of negli- 
gible shear and in a layer whose Mach number is almost constant and is 
less than that produced by a sharp cone or wedge. 

Of particular significance is the fact that the region of reduced 
Mach number near the surface is also a region of reduced Reynolds number.^ 
Until the boundary layer engulfs this region, its stability and transi- 
tion characteristics, as well as its flriction and heat-transfer charac- 
teristics, should be those associated with the reduced Reynolds number. 
This reduction in Reynolds number near the surface of blunted bodies 
explains the downstream movement of transition observed in references 1 
to 3, and is the basis used in this report for comparing the boundary- 
layer characteristics of blunted and unblunted bodies. 


■^Tiis reduction in surface Reynolds niunber due to blunting and its 
effect on laminar heating have recently been independently calculated in 
ref. 6 for hypersonic speeds. Ro attempt was made, however, to define 
the thickness and axial extent of the low Reynolds number layer or its 
effect on transition location. 
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AMLYSIS 

The Mach number in the Inviscid shear layer produced near the sur- 
face of blunted cones and wedges Increases continuously from the surface 
value to the value that would exist at the surface of the corresponding 
unblunted bodies. The Reynolds number per unit length at the outer edge 
of the boundary layer ("outer-edge" Reynolds number) therefore remains 
less than the free-stream (or unblunted) value until the boundary layer 
absorbs the entire shock-produced shear layer. If the transition point 
is determined primarily by the Reynolds number at the outer edge of the 
boundary layer, a progressive downstream movement of transition would 
therefore be expected as the leading edge or tip bluntness is gradually 
increased. The maximum downstream movement would be expected when the 
blunting is sufficiently great so that the outer-edge Reynolds number 
is close to the Inviscid surface value for the entire laminar run. In 
the following sections, the maximum reduction in outer-edge Reynolds 
number is calculated, and the blunted area required to produce this 
maximum reduction over the entire laminar layer is estimated. 


Reduction in Surface Reynolds Rumber Due to Blunting 

At a station sufficiently far downstream of the vertex, where the 
surface static pressure for a blunted body is close to that of the un- 
blunted body, the Reynolds number near the surface can be written as 

_ Pn^n /■ 

Rei = Ml 


where subscripts n and 1 refer to inviscid surface values for the 
blunted and unblunted bodies, respectively. (All symbols are defined 
in appendix A.) These inviscid surface values will be assumed, as 
usual, to represent the outer-edge conditions that determine boundary- 
layer development. 


The use of Sutherland's viscosity equation yields 

.2 


Re^ 


_ AlV /^ + 

" W/ W + S/% 


( 2 ) 


Dividing the numerator and denominator by the ambient static tesmperature 
tQ and converting to Mach number functions yield 






( 3 ) 
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( 5 ) 


The inviscid surface Mach numher for the "blunted body is de- 
termined by the ratio where is the static pressure at the 

surface of the unblunted body and is the stagnation pressure down- 

stream of a normal shock at the free-stream Mach number Mq. The in- 
viscid surface Mach nxunbers are shown in figure 1 as a function of Mq 
for several cone and wedge angles. Since the total pressure is 

less than the total pressure at the surface for unblunted bodies, the • 
surface Mach number for the blunted bodies is less than the sur- 

face Mach number for the unblunted bodies M^. The difference between 
Mq and M^^ increases as flight Mach number increases . 

The Reynolds number ratio of equation (3) is plotted in figure 2 
for the same cone and wedge angles as those in figure 1. This ratio 
decreases rapidly as flight Mach n\miber increases. If the transition 
Reynolds number is unchanged when the leading edge or tip is blunted, 
and if the blunting is adequate to coyer the laminar boundary layer with 
a sufficiently thick layer of low Rejraolds niunber air, then it should be 
possible to increase the length of laminar run by a factor inversely 
proportional to the Reynolds niimber ratio of figure 2. For slender 
cones and -vfedges, the possible increases in laminar run range from fac- 
tors of the order of 2.0 at Mq = 3.0 to 10 at Mq « 8.0 and 30 at 

Mq =< 15.0. The significance of such large increases in laminar run for 
reducing the heat-transfer rate and friction drag for very high-speed 
aircraft is apparent. 

Evidence that increases in laminar run of the magnitude indicated 
by figure 2 are actually attainable is presented in references 1 to 3. 
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In reference 3, for example, the transition point at Mq = 3.1 -was 
moved downstream by a factor of 2 (from 5 to 10 in. at a Reynolds 

number of 3.56xi0^/ln .) when the leading-edge thickness was increased 
from 0.0008 to 0.008 inch. This experimental movement of transition 
compares very favorably with the VEilue 2.17 predicted on the basis of 
the Reynolds number reduction shown in figure 2(b). 

In reference 2, downstream movements by factors ranging from 2.3 
to 3.6 were observed for a blunted flat plate at various angles of 
attack at Mq = 4.04. The movement predicted by figure 2(b) for this 

Mach number is 3.57. For swept wings, little or no downstream movement 
was observed in reference 2. This is in agreement with the expected 
weakening of the leading-edge shock due to sweepback. Whether down- 
stream movements of the order of magnitude predicted by figure 2 are 
attainable at higher Mach mnaber or for other body shapes remains to 
be established by further experiments. 


Estimation of Bluntness Required to Obtain 
Maximum Movement of Transition 

In order to determine the bluntness area required to cover the 
entire laminar boundary layer with a low Reynolds number layer of neg- 
ligible gradient, it is convenient to define a thickness of this layer 
which limits the Mach nianber to veilues near the inviscid surface value. 
A suitable thickness is the distance from the surface to the streamline 
that passes through the sonic point of the detached shock wave (point 
where the Mach number Just behind the shock is unity) . From the vertex 
to the sonic point the stagnation pressure downstream of the shock does 
not vary greatly; consequently, the Mach number should remain near the 
inviscid surface value in the layer thus defined. 

An expression for the thickness iriJJL be derived under the assump- 
tion that the shear profile produced by the detsiched shock does not 
diffuse or dissipate, that is, the profile remains unchanged until it 
is engulfed by the boundary layer. The rate of dissipation of the 
shock -produced shear layer is discussed in appendix B, and tends to in- 
crease the bluntness area required to produce a given thickness of the 
low Mach number layer. 

With dissipation neglected, the thickness of the low Mach number 
region can be estimated by means of the detached-shock-wave theory of 
reference 7. In this theory, the detached shock wave is assumed to 
have a hyperbolic form independent of the shape of the body that pro- 
duces it. This form has been found to a^ee well with experimental re- 
sults for a large range of body shapes in the moderate supersonic Mach 
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number range (ref. 8), but becomes more questionable as the Mach number 
increases. The portion of the shock from the vertex to the sonic pointy 
however, can be satisfactorily represented by the assumed hyperbola for 
all Mach numbers. This fonn of the shock wave should, therefore, yield 
a satisfactory estimate of the thickness of the low Mach ntnnber layer as 
well as the shock location, to the extent that the other assumptions of 
the theory (constant specific heat, inviscid flow, etc.) are valid. 

In order to estimate this thickness, the Mach number in the layer 
is assumed to be constant at a value corresponding to the mass centroid 
of the layer. This Mach number, denoted by Mq, is determined from the 
ratio p^/Pq^ where p^ is the static pressure on the surface of the un- 
blunted cone or wedge and Pq is the total pressure downstream of the 

shock on the centroid streamline- (A simple and satisfactory estimate 
of Pq can be obtained by using the arithmetic mean of the stagnation 

pressures at the sonic point and at the vertex.) The continuity equa- 
tion for the layer shown in sketch 2 can be written 

PcA^(aVa)m(, - PoAs(a*/a)mo (6) 



where A^ is tbe area of the low Mach number layer and Ag is the free- 

stream area of the stream tube between the axis and the shock sonic point 
S. If the bluntness of the body is defined as its cross-sectional area 
at the sonic point Agp (sketch 2), the ratio of the area of the low 
Mach number layer to the blunted area becomes 

A _ As Pq 

%B Agg Pq (aVa)]^ 


( 7 ) 
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From reference 1 , 


As 1 

Agg X — B cos T] 


where ® ^ (A*/a)mq, and r\ 

line defined in reference 7 . 


is the mean inclination of the sonic 


( 8 ) 


The area of the low Mach number layer defined by equation (7) is 
shown in figure 3. This area is seen to increase rapidly with increas- 
ing Mach numbers for the blunted flat plate (8y = O) and for the blunt- 

nosed cylinder (8^ = O) . However, for wedge half -angles greater than 

5° and cone half-angles greater than 10°, the area does not vary greatly 
with Mach number. 


With the thickness of the low Mach number and low Reynolds number 
layer thus defined, the blunting required to provide a low external- 
stream Reynolds number for the entire laminar boundary layer to the ex- 
pected or hoped-for transition point can be estimated. This is done by 
calculating the laminar boundary-layer thickness at the expected tran- 
sition Reynolds number, which is based on conditions in the low Mach 
number layer near the siarface. By equating this thickness to the thick- 
ness of the low Reynolds number layer, the required values of the blunt- 
ness area can be calculated. 


Thus, for blunted wedges the required ordinate at the body sonic 
point is 

^tr 

wMle for blunted cones (with Str« ^l)^ 


( 9 ) 


ysB 


■,tr^tr 

j/Asb) J 


1/2 


( 10 ) 


where A^/Agg is given In figure 3, and r-j^ is the radius of the 

blunted cone at the expected transition point. Equations (9) and (lO) 
show that the amount of blunting required to cover the laminar boundary 
layer with a low Reynolds number layrar is not large. For the wedge 
(eq. ( 9 )), the ordinate of the body at the sonic point need be only of 
the order of magnitude of the boundary- layer thickness at the expected 
transition point} for the cone (eq. (lO)), the required radius of the 
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body at the sonic point is of the order of the geometric mean of the 
body radius and the boundary- layer thickness at the expected transition 
point. 


The required bluntness is considerably reduced if the displacement 
effect is considered, since the low Reynolds number layer is moved away 
from the surface by an amount equal to the displacement thickness of the 
boundary layer (ref. 3) . The required values of ygg should therefore 

be calculated with (6 - 5*)^ in place of &tr in equations (9) and 

(lO) . Expressions for 5 and &* for constant surface temperature 
were obtained from equations (18) and (22) of reference 9, based on the 
flat-plate theory of reference 10. The value of 5 was assumed to cor- 
respond to —a 0-99. At the transition point, these expressions ’can 

Uoo _ — - 

be combined to yield (for y = 1.40) 



( 11 ) 


where Ltr Is the distance along the su rface to the transition point 
and C is the proportionality constant in the linear viscosity- 
t^perature variation. For cones, this expression is divided by -y/3. 
In terms of (S - 6*). , the~bluntness required to cover the laminar 

wl* 

boundary layer with a low Reynolds number layer becomes 


ySB (^A*)tr ■ (^*A)tr 

^r 


for wedges and 


Ltr 




1 

tr 

fs). ■ ' 

L tr 



1 (-^Asb) j 


1/2 


( 12 ) 


(13) 


for cones . 

With equation (l2), the calculated bluntness areas agree as closely 
as could be expected with the experimental values that produced the 
maximum downstream movement of transition in the e:q)eriments of refer- 
ence 3. This maxim um downstream move ment was found jbp^tgke pleLce,.for a 
leading-edge thickness of about 0.008 inch, which is about two-thirds 
of tE^ calculated value . Further 'indreases in leading-edge thickness 
had no appreciable effect on transition location. Since the thickness 



10 


MCA TN 3653 


of the low Reynolds number layer has heen rather arbitrarily defined^ 
the values of ygg calculated from equations (12) and (l3) can he re- 
garded only as rough estimates of the blunting required to produce 
maximum downstream movement of transition. These values should, how- 
ever, be on the conservative side, since they produce nearly the maixi- 
mum possible reduction in Reynolds number and Mach number over the en- 
tire laminar boundary layer. 


Invlscld Mach Number Profiles for Blunted Cones and Wedges 

Although the maximum effect of blunting on boundary- layer develop- 
ment and treLnsition depends on the portion of the shock-produced shear 
layer near the surface, the entire shear profile is of interest if the 
outer edge of the boundary layer moves out of the low Reynolds number 
layer defined In the preceding section. In order to determine the na- 
ture of the entire shock-produced shear profile, the shape and location 
of the shock must be prescribed. For moderate supersonic speeds, the 
hyperbolic foim assumed in reference 7 is adequate; but as the flight 
speed approaches the hypersonic range, the portion of the shock beyond 
the sonic point is increasingly Influenced by body shape. This situa- 
tion arises partly because the region between the shock and the body 
beccanes smaller as Mq increases; consequently, characteristics from 

portions of the body far downstream of the sonic point reach the shock 
before it has decayed to its asymptotic strength. .In addition, the 
overexpansion near the shoulder of slender blunted bodies , which takes 
place at lower speeds, gradually becomes an underexpansion at hypersonic 
speeds, that is, a Prandtl-Meyer expansion from the sonic point fails 
to reduce the pressure to, or below, the asymptotic static pressure. A 
rather long process of reflection of expansion waves between the shock 
and the body must, therefore, take place before the asymptotic pressure 
is reached on blunted cones or wedges. 

This consideration also affects the distance required to obtain 
the inviscld surface Mach numbers and Reynolds numbers calculated in 
the preceding sections. A more accurate evaluation of the effect of 
blTintlng would include the variation of outer-edge Mach number and 
Reynolds number along the entire body due to the pressure gradient. 
Perhaps mean values of these nxmibers could be used to predict the loca- 
tion of treinsitlon. These mean values would be lower than those shown 
in figures 1 and 2, which means that the predicted transition point 
would be even farther downstream than if the pressure gradient is neg- 
lected. The fact that the self-induced pressure gradient is entirely 
favorable also tends to increase the stability of the laminar layer. 

The effect on transition of increasing flight Mach number thus appears 
to be a favorable one. For computing the shock-produced shear profile, 
however, the narrowing region between shock and body as flight speed 
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increases introduces difficulties, in that no general shock shape is 
available beyond the sonic point, and the asymptotic profile may be so 
far downstream as to have no practical significance. It was neverthe- 
less felt to be worthwhile to compute these asymptotic profiles for very 
high Mq with the hyperbolic shock form of reference 7, if only for 
comparison with more accurate future computations based on experimented 
shock forms or exact characteristic solutions for particular bodies . 

The variation with Mach number of the hyperbolic shock of reference 7 
is similar to that which would be expected, in that it decays niuch more 
slowly for hypersonic speeds than for moderate supersonic speeds. 

The computation method is presented in appendix C, and the result- 
ing asyr^jtotic invlscid shear profiles are shown in figure 4 for flight 
Mach numbers from 2 to 20. Indicated on each profile is the thickness 
of the low Mach number layer as defined in the preceding section. It 
is seen that this definition does, in fact, restrict the Mach number to 
values close to the surface value. 

The profiles for blunted wedges differ qualitatively from those of 
blunted cones at all Mach numbers. For the blunted wedges, the Mach 
number gradient is zero at the surface; whereas, for blunted cones the 
gradient has a positive value. Since these gradients depend on the form 
of the shock near its vertex, they should be correct for all Mach num- 
bers ^■rithln the limitations of the other assumptions of the analysis 
(static pressure equal to values for the unblunted body, constant spe- 
cific heat, inviscid flow, etc.). The portions of the profiles above 
the boundary of the low Mach number layer should be good approximations 
for Mq less than 5.0, but seem to become much too thick at higher Mach 
numbers, particularly for the flat plate (e^ = O) and for the blunt- 
nosed cylinder (d,, ~ O) . This thickness is associated with the very 
slovr decay rate of the assumed hyperbolic shock at these Mach numbers. 
Since the shock lies quite close to the body at hypersonic speeds, these 
profiles would, as previously surmised, be applicable, if at all, only 
at very large distances from the vertex, where the shear layer is thin 
compared with the distance from the surface to the shock wave. At such 
distances, of course, the boundary layer, which is of the same order of 
thickness as the layer between the shock and the body at hypersonic 
speeds, would already have engulfed the entire shock-produced shear 
layer . 

For higher cone and wedge angles, the shock decays more rapidly to 
the asymptotic strength, and the resulting profiles appear to be more 
in harmony with expectations . 

Although the conrputed profiles beyond the boundary of the low 
Reynolds number layer are not reliable at high Mach numbers, they agree 
well with measured profiles at Mach 3.1 (ref. 3). If more accurate 
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shock-produced profiles are desired for higher speeds ^ the shock fonn 
must he calculated for eeich body shape. Such computations vould be use- 
ful for estimating the variation of transition location as the blunted 
area is gradually increased, but are not required for estimating the 
maximum downstream movement, or the blunted area required to produce 
this movement . 


Changes in Laminar Recovery Temperature, Heat-Transfer Rate, 

and Friction Coefficient 

The downstream movement of transition due to blunting means that 
larger portions of the aircraft surfaces will be subjected to laminar, 
rather than turbulent, heat-transfer rates and friction coefficients. 

The blunting should, therefore, produce substantial reductions in over- 
all heat-transfer rate and friction drag. There is, however, an in- 
crease in laminar equilibrium recovery temperature corresponding to the 
reduction in Mach number, and a change in laminar heat ■^transfer rate and 
friction coefficient due to the reduction in Reynolds number. These 
must be evaluated in order to estimate the magnitude of the advantages 
due to blunting. 

The heat-transfer coefficient and friction-drag equations of refer- 
ence 11 are used for this estimate. Although these equations are based 
on the assumptions of constant specific heat and Prandtl number , and no 
dissociation, they agree in order of magnitude with more exact numerical 
computations even at hypersonic speeds (ref. 12). 

The ratio of laminar heat-transfer rate with and without blunting 
can be written 


% 

<ll 




(14) 


where D and r ai’e defined by equations (4) and (5), tg^n eind te,l 
are equilibrium recovery ten5)eratures with and without blunting, and 





is the shear function of reference 11 evaluated at the sur- 


face. This function is given in reference 11 for several outer-edge 
Mach numbers Mo, and several ratios of surface temperature to outer- 


edge temperature 


t^/teo* 
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The temperature-difference ratio in equation (14) can be written 


'*'e,n " ■^w “^0 

1 + 0.845| 

r 3 

. 1 + 0.2M^ 

w j 

- . 2 1 
^0 1 + 0.21<^ 

^1 + 0.169M^^ 

1 - S 

1 to 

■^0 

B 

§ 

1 + 0.845| 

l' ^ 

_ 1 + o.2M^ 

^0 2 * 
1 + 0.2M^ 

^1 + 0.1691^^ 

■*=w 
1 - t^ 
^0 

(15) 


where the laminar recovery factor is assumed to he 0.845 and y is 1.4. 

The ratio of equilibrium surface temperat\u:*es with and without 
blunting is obtained from the definition of recovery factor: 


whence 


^e,n “ ^ ^e,l ^1 
To - tn “ To - t 



0.155 
1 + 0.2J^ 
0.155 
1 + 0 . 21 ^ 


0.845 


0.845 


(16) 


(17) 


The laminar-skin-frlction ratio is, firom the equations of reference 

11. 



The ratios of laminar recovery temperature, skin friction, and heat 
transfer for flat plates are shown in figure 5. Although there is a 
slight increase in laminar equilibrium temperature for the blunted flat 
plate (this was observed experimentally in ref. 3) , the laminar skin 
friction is reduced over the entire range of flight Mach ntnnbers , and 
the heat-transfer rate is reduced except for wall temperature near 
equilibriiim. (The rapid increase in the heat-transfer ratio near re- 
covery teii5)erature arises from the small increase in recovery temperature 
due to blunting. The heat transfer without blunting approaches zero 
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for these values, whereas the heat transfer with blunting becomes smal 1 
but is not yet zero.) Figure 5 shows that blunting the leading edge of 
a flat plate or cylinder can produce, in addition to the longer laminar 
3Tun, a small but significant reduction in the skin friction and heat- 
transfer rate of the laminar boundary layer itself. 


Cooling Requirements for Stability 

The static temperature at the edge of the boundajTr layer is con- 
siderably higher for a blunted cone or wedge than for sharp bodies. 

The ratio tu/tp is, in fart, given by i/d (eq. (4)). This increase 
in outer-edge temperature means that, for a given surface temperature 
tv, the ratio tv/tn is smaller than tv/ti. The outer-edge Mach num- 
ber is also reduced. Shown in figure 6 are the outer-edge conditions 
for a blunted and unblunted flat plate, and for a blunted and unblunted 
10° half -angle cone for a surface-to-ambient temperature ratio of 4.0. 
These conditions are compared with two of the laminar stability limits 
given in references 12 and 13. This comparison shows that blunting 
moves the outer-edge conditions far into the stable region in the hyper- 
sonic speed range. (Although the stability-range curves shown are based 
on two-dimensional disturbance theory, recent conputations by Dunn and 
Lin (ref. 13) indicate that three-dimensional disturbance theory also 
yields laminar stability to extremely hi^ Reynolds number but that some- 
what lower surface temperatures are required.) 


Effect of Blunting on Heat-Transfer Rate Near the Nose 

In order to estimate more accurately the net decrease in heat- 
transfer rate due to blunting, it is necessary to determine how the 
heat-transfer rate near the nose of the blunted cone differs from that 
on the pointed cone. An estimate of this difference can be made by 
comparing the heat-transfer rate for the sharp conical nose with that 
for the inscribed spherical nose (sketch 3) . 



Sketch 3 
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This heat-transfer ratio can he -written as 

%p _ ^sp ^sp (*^0 ~ 

Ac (Te^i - t^) 


(19) 


■where the subscripts sp and c refer to the spherical and conical 
nose; respectively. .The area ratio of equation (19) is 


V ^ ^c 

Ac cos 0c 


(1 - sin 6c) 


( 20 ) 


and the tempera-ture-difference ratio is 


*^0 - “^w 
^e,l “ ^ 


1 + 


r - 1 




^=0 


1 + 


r - 1 


“o 




° / 

IV 


r - 1 





( 21 ) 


The mean heat--transfer coefficient for the spherical nose is asstuned to 
be the stagnation -point value presented in reference 14: 



( 22 ) 


where c = — a / - — C g.^ and (Nu/-y/Re^) is about 0.61 for 


^ / f2_ 

r ^ '"P;.St W'sp 

Prandtl mmber of 0.72 and for a ratio of wall tempera-ture to stagnation 
temperature (t^/To) close to zero (corresponding to cooled surfaces at 


very high Moj.^ The stagnation pressure coefficient 1-84 for 

y = 1.4. The mean cone heat-transfer coefficient is ; from reference (14); 



(23) 


2 

Since the publication of ref. 14, Reshotko and Cohen have found 
that the expression for c given therein for supersonic flow is in 
error. The correct expression for this constant is that given above. 
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where (Nu/-y/5e^^ is 0.51. The ratio of heat-transfer coefficients, 
therefore, hecomes 



where Pg^ is the stagnation pressure behind a normal shock at Mach 
number Mq. The heat transfer ratio obtained by substituting equations 

(20), (2l), and (24) into equation (l9) is plotted in figure 7 for cone 
half-angles of lO^^and 20° and for a surface-to-ambient static- 
tenperature ratio of 1.0. The over-all heat transfer for the inscribed 
spherical nose is seen to.be less than half as great as that for the 
conical nose. The blunted nose, therefore, has the advantage of a lower 
heat- transfer rate near the vertex as well as along the downstream 
surfaces. 


DISCUSSION 

The preceding sections have shown that the Reynolds number per unit 
length at the outer edge of the boundary layer is lower for blunted fuse- 
lages and wings than for unblunted ones. The limited data available 
agree with the conclusion that the transition location can be increased 
by a factor of the order of -the ratio of the surface Reynolds number 
without blunting to the surface Reynolds number with blunting. This 
factor increases rapidly with increased flight speed , particularly for 
moderately slender wings and bodies. 

As- an example of the magnitude of this effect, a 10° half-angle 
cone at a Mach number of 15 -trill be considered. If the transition point 
is located 1 foot do-wnstream of the vertex without blunting, it might, 
on the basis of figure 2, be moved 25 feet downstream of the vertex if 
the tip is blunted. 

The bluntness required is, from equations (U) and (13), 

^ . o.sU-Y''* 

^ \^®tr/ 

If ( - - - ) is of the order of 10"^, then the required value for ygg is 

^°try 

only about 5 inches. The ratio of the blunted area to the cross-sectional 
area of the cone at the transition point is, therefore, approximately 
0.01. If the -transition point (25 ft) is near the end of the body, the 
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over-all heat-transfer rate would he reduced hy blunting from the value 
corresponding to almost completely turbulent flow to the value corre- 
sponding to completely laminar flow. The blunted cone would, therefore,, 
heat up much more slowly than the pointed cone and would require much 
less coolant to maintain a given surface tenperature. The ratio L^.j-/yn 

is about 300^ therefore, the effect of dissipation of the shock-produced 
shear layer can probably be neglected (see appendix B) . 

Furthennore, during the heating process the ratio of surface tem- 
perature to outer-edge ten 5 )erature remains much lower for the blunted 
cone or wedge (fig. 6) so that the adveintages of cooled surfaces with 
regard to laminar stability prevail longer than for the pointed cone or 
sharp-edged wedges. Both the low surface Eeynolds number and the higher 
outer-edge temperatures work toward preseirvation of laminar flow for a 
much larger distance along the. surfaces of blunted bodies and wings. 

These advantages with regard to increased laminar run and Increased 
laminar stability appear to involve no serious disadvantages. The fric- 
tion drag is reduced, and the total drag should not increase appreciably 
for the small required values of the bluntness ratio. Reference 16 shows 
that, for spherical -tipped cones of fixed total length, the total drag 
to Mach number 7.0 is very near the value obtained for the sharp-tipped 
cone for ratios of nose diameter to maximum body diameter less than 0.25. 

The quantitative effects of blunting on transition location pre- 
viously computed are based on the hypothesis that the transition Reyn- 
olds mnnber is substantially unchanged when a body with a sharp tip is 
blunted. Although this hypothesis produces good agreement with the 
experimental results of reference 3, the possibility should cejrtainly 
be kept in mind that, at higher Mach numbers or with other body shapes, 
the transition Reynolds mmiber may be altered by such factors as pres- 
sure gradient and outer-edge Mach nijmber and Reynolds number. Further- 
more, as the length of laminar run increases, the possibility of pre- 
mature transition due to surface roughness or stream turbulence also 
increases, and the dissipation of the shock-produced shear profile be- 
comes important . Whether any of these factors will seriously reduce the 
attainable downstream movement of trsinsitlon due to blunting remains to 
be determined experimentally. 

Many theoretical, problems also require solution before the quanti- 
tative effects of blunting on transition can be predicted with confid- 
ence. One basic problem, of course, is that of the development of a 
laminar boundary layer in a nonuniform external stream. Solution of 
this problem would establish the magnitude of external shear that is 
negligible and, consequently, the conditions for which the boundary 
layer can be assmed to develop in a layer of reduced Reynolds number 
corresponding to the mean value near the surface. This solution might 
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reveal whether, as indicated by the results of reference 3, it is suf- 
ficient, in general, to obtain the maximum reduction in Eeynolds number 
only for the inner half or trwo-thlrds of the boundary layer rather than 
at the outer edge. The latter question, however. Involves predicting 
the location of transition for various velocity profiles, which cannot 
as yet be done even for laminar layers in a uniform external stream. 
Since the required blunting is small, however, this question appears to 
be of secondary importance. 


CONCLUDING EEMAEKS 

It is clear from the preceding discussion that many questions re- 
main unanswered in this report. The principal observation that the 
Reynolds nvunber and Mach number near the surface are reduced by blunt- 
ing, and also the approximate meignitude of the reductions are fairly 
well established. The assumption that the boundary-layer development 
should be determined primarily by the reduced Eeynolds number and Mach 
niimber near the surface rather than by the flo\^ outside the inviscid 
shear layer also seems reasonable. The principal benefits from blunt- 
ing, however, lie in the hypersonic speed range, where many of the quan- 
titative results calculated herein are subject to corrections whose mag- 
nitude is as yet unknown. Qualitative estimates indicate that some of 
these corrections, such as the displacement effect or the pressure gra- 
dients, either inviscid or self -induced by the boundary layer, should 
have a favorable effect on the downstream movement of transition. Other 
effects, such as surface roughness, stream turbulence, or changes in 
transition Reynolds number, may tend to limit the downstream movement 
of transition to values less than those predicted. Dissociation at 
very high Mach numbers may have a significant effect on outer-edge con- 
ditions and, consequently, on the maximum transition movement to be ex- 
pected. As usual, when so many unknown factors contribute to a phenom- 
enon, experiment must be relied upon to determine which factors are 
dominant and which are of minor importance. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, November 21, 1955 
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APPENDIX A 
SYMBOLS 


The following symbols are used in this report: 
area 


(A* /a)j^ 


B 

C 

C 


PjSt 


isentropic area contraction ratio from Mach number M to 
Mach number 1.0 

speed of sound 


I (A»/A) 


Mo 

constant in linear viscosity-temperature relation 
stagnation pressure coefficient 


U, 


_0 

r 


PQ 

Pst 


Cp^st 


D 



h 

k 

L 

M 

Nu 

P 


1 + 0.2M^ _ t3_ 

1 + 0.2J^ "^n ■ 

shear function at surface (ref. ll) 
heat-transfer coefficient 
thermal conductivity of air 
length of conical tip 
Mach number 
NuBselt number 
stagnation pressure 


P 


static pressure 
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heat-transfer rate 
Reynolds numher 
recovery factor or radius 

cone radius at station where profiles are determined 

Sutherland's constant for air, 198.6° R 

stagnation temperature 

static teng)erature 

surface temperature 

velocity 

velocity downstream of normal shock ahead of spherical 
nose 

distance along surface 
coordinate normal to surface 

'\|“S - ^ 

tQ + S 

t]_ + s 

ratio of specific heats, 1.40 
boundary- layer thickness 
boundary- layer displacement thickness 
inclination of sonic line (ref. 7) 
semivertex angle of cone 
semivertex angle of wedge 


coefficient of viscosity 
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P 

T 

fp' 


density 

shear force at surface 
shock angle 


Subscripts : 


C 

c 

e 

n 

S 

SB 

sp 

st 

tr 

W 

V 

0 

1 


centroid 

cone 

equilibrium 

inviscid surface values for blunted cones or wedges 

sonic point oh detached shock wave 

sonic point on body 

sphere 

stagnation 

transition point 

wedge 

surface values 
ambient conditions 

inviscid surface values for unblunted cones or wedges 
value at outer edge of boundary layer 


Superscript : 


1 


local conditions in inviscid shear layer 
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APPENDIX B 


DISSIPATION OF THE LOW MACH NUMBER LAYER 

In order to estimate more closely the amount of blunting required 
to maintain a given thickness of the low Mach number layer, the rate of 
dissipation of the shock -produced profile must he considered. The sim- 
plest method for estimating the rate is to consider the profile produced 
by the detached shock wave as a step function (sketch 4) : 



in which the outer velocity is that corresponding to the unblunted body 
and the inner velocity is that produced near the surface by blunting the 
vertex. The profile dissipation can then be considered identical to 
that at the Interface of two parallel laminar jets emerging at the same 
static pressure. The equation for the velocity profile in the inter- 
action region is given in reference 17 for the case when Uj_ and 
differ by a small amount . The appropriate equation is : 



where y^ is the initial thickness of the lorii Mach ntmiher layer , ^^®y^ 

is Reynolds number based on y^ and outer-flovr conditions, and i>(a) 
is the enror function of a. Profiles calculated from equation (B1) for 
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^ r 

Rey^ = 10 and for several values of x/yn are shown in figure 8. 

Apparently, the velocity near the surface does not change appreciably 
until x/yjj is of the order of 1000. Although these profiles are 

valid only for small differences between U 2 _ and the order of 

magnitude of the dissipation remains the same for large differences 
(see fig. 4.11, ref. 17). The value of x/y„ < 1000 is, therefore, 

probably a good estimate for the length of run in which dissipation 
of the shock-produced shear profile can be neglected if this profile 
remains laminar . 


If transition to turbulence takes place in this layer, the length 
of run for which dissipation can be neglected is appreciably reduced. 

No experimental results are available to estimate under what conditions 
the shock-produced shear profile is likely to undergo transition. How- 
ever, an indication of whether transition is a possibility in this 
layer can be obtained from the stability criterion for parallel jets 
developed in reference 18. This criterion states that the interface 


can become unstable if the quantity 




P 


vanishes in the interface 


profile. However, the profile is stable if this quantity vanishes only 
at points in the profile where the velocity satisfies one or both of 
the following Inequalities : 


U < - a^ 

U > + a„ 


These conditions sissure that disturbances ffcm either stream will not 
reach a layer in which amplification is possible. 

In terms of Mach number profiles, these conditions can be stated 
as follows : 


The profile is stable to two-dimensional disturbances if the quantity 


_d 


dM/dy 




1 + ^ I 


vanishes only at points where 


Ml 


M < 


V(r - i)(Mi + 2) 


(B2) 
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or where 


Mn + 1 

“ (r - 1)(2 - M^) 

The latter condition cannot he satisfied for > 2.0. 
the stability of the profile depends chiefly on whether 


vanishes only where condition (B2) is fulfilled. Some sample computa- 
tions based on the profiles of figure 4 indicate that condition (B2) is 
generally satisfied for the hlunted-cone profile hut not for the hlunted- 
wedge profiles. The latter profiles therefore are more inclined to 
undergo transition than the former. If transition occurs^ the amount of 
hluntness required to produce a prescribed thickness of the low Mach num- 
ber layer at a given station may be considerably greater than calculated 
on the basis of laminar flow. This discussion must necessarily be incon- 
clusive, since the location or even the existence of transition cannot be 
established from stability theory alone. 


(B3) 


Consequently, 
d P dM/dy 

<3y 


1 + 
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APPENDIX C 


MACH HUMBEE PROFILES PRODUCED BY DETACHED SHOCK WAVES 

Tb.e Mach, number profiles normal to the surfaces of cones and wedges 
with slightly blunted tips or leading edges can be calculated from the 
one-dimensional continuity equation if (l) the form and location of the 
detached shock wave is kno™, (2) the static pressure is constant normal 
to the surface^ and (3) diffusion and dissipation of the profile are 
neglected. Condition (l) is most conveniently satisfied by using the 
detached-shock-wave theory of reference 7. Condition (2) is satisfied 
at stations sufficiently far downstream of the nose or leading edge, 
where the surface static pressure has reached, or closely approached, 
the value obtained with unblunted cones or wedges. At moderate super- 
sonic Mach numbers, the required distance is of the order of 3 to 10 
times the thickness of the blunted portion of the nose or leading edge. 
This condition is not quite satisfied for blunted cones, because the 
flow field approaches a conical distribution characterized by a gradual 
decrease of static pressure from the surface to the shock wave. But if 
the profile extends only a small portion of the distance from the sur- 
face to the shock wave, this gradient can be neglected without serious 
error. Condition (3) remains an assumption whose validity decreases as 
the distance along the body increases. It implies that the profile re- 
mains unchanged in form for an unlimited distance downstream of the ver- 
tex. As pointed out in appendix B, this assung)tion appears to be fairly 
good for distances of the order of 1000 times the thickness of the 
blunted portion of the body if the profile remains laminar. 

The profile computation is set up with the aid of sketch 5, .which 



Sketch 5 
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applies either for cones or wedges. If the sonic-point area of the body 
is used for reference purposes, the continuity equation can be written as 


or 




(a* /^)Mq 


P* (a*/A)m' 



(Cl) 


(C2) 


where the primes refer to local conditions in the inviscid shear layer. 
Since the stagnation pressvire along each streamline remains constant 
downstream of the shock, P'/Pq is the stagnation-pressure ratio across 
the shock at the point where the streamline bounding the area A enters 
the shock. If the shock angle at this point is (p, then the total- 
pressure ratio can be written (ref. 19) 


P' 


P 


0 


6M^sin^<p 
, M^sin^<p + 5 


7/2 



(C3) 


The Mach number M' at the area A* can also be expressed in terms of 
shock angle by the relation 


1 


+ 0.2M'^ =■ 


/pA ^ 

W) ~ V^o Po Pi/ 


(C4) 


The function (A*/a)jji, as well as P'/Pq, is a function of the shock 

• angle «p at the point where the streamline crosses the shock wave. The 
differential d(A/Agg) of equation (C2) must be converted into a func- 
tion of <P in order that the integration may be carried out from 
<P = 90^ to cp = qq, where qq is the shock angle corresponding to the 
unblunted cone or wedge. The Mach niimber M' as a function of A' can 
then be obtained from equations (C4) and (C3). 
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From reference 7 (eq. (5)), the relation between shock angle and 
shock ordinate for the assiimed hyperbolic wave is 


where 



(xo/ysB) 


2 


p2(p2^an2(p - l) 


fo_ 

^SB 


3 ^r~ a / P^tan^cpg - 1 

^SB V ^ 


(C5) 


(C6) 


and (pg is the shock angle at the shock sonic point. The ratio yg/ygE 
is a function of Mq and depends on whether the flow is two-dimensional 
or axially symmetric. 

The area differential of equation (C2) can now be expressed as 
follows : 


For two-dimensional flow: 


A \ ^ y \ _ Q fo_ tan q> sec y d(p 

\^By [yssj ysB 


For axially symmetric flow: 


/ /y 

\2 1 

~ p " 



tan q) sec (p dtp 

VSB J “\ySB/ 

(p^tan^T - l)^ 


(C7) 


(C8) 


Combination of equations (C2) to (C8) yields the following final 
expressions for determining the variation of M' with A' : 


For plane flow: 


A’ _ y' 
Asb ysB 






1/2 

Po 

A + 0.2Mg > 

1 — f- 

Pi 

\1 + 0.2M’^ 

1 M- 

V. 

90° 



tan 


(p^tan^ 




(C9) 
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For axially-symmetric flow: 


( 2 2 \P0 P fl + 0. «0/ tag <t> sec^<P \ 

VSB/ \ ® jpi / \1 + 0.2M'^ ^ \(p2tan^ -1)2/ 

^0° (n 


d<P 


A' ■■ „2 

(CIO) 

In these equations, y' is the linear distance normal to the wedge or 
cone. For y’ « rq, the area ratio A>/Agg in equation (lO) is equal 
2r]_y* 

to -2 

^SB 


Equations (C9) and (CIO) have heen integrated numerically for sev- 
eral Maeh numbers and for several wedge and cone angles. The resulting 
Mach number profiles are shown in figure 4. 
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Figure 1. - InviBcid surface Macli number for 'blunted cones and tredges. 
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Fll^xt Macli number^ Mq 

(a) Cones. (t) Vei^s 

Flgore 5. - Area of lov Mach nnaiber layer for hlimted cones and vedges. 
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1.2 1,6 2,0 
Mach number^ M 


(h) FU.^t Mach nuniber, 2,0j cones. 

Figure 4, - Irnriseid Mach number prof lies for hlunted cones and vedges. 
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(r) •p'H gVit . Mach number, 3.0; vedges. 
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r-l 

U 



.6 1.0 1.4= 1.8 2.2 

Mach number, M 

(d) Fll^t Mach number, 3.0j cones. 

Figure 4. - Continued. Inviscid Mach number profiles for blunted cones 6ind 
vedges. 
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(e) FUgJit Mach number, 5.0j Vedges. 
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2 3 

Mach number, M 

(f) Flijfi't Mach number, 5.0; cones. 

Figure 4. - Continued. Inviscid Mach number profiles for blunted cones and 
vedges . 
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200 




100 



Fh 


(n) Fli^t Mach number, 20.0; cones. 

Figure 4. - Concluded. Invlscid Mach number profiles for blunted 
cones and wedges . 
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Outer-edge Mach number, 


(a) 10° Half -angle cone. 


(h) Plat plate. 


Flpre 6. - Effect of blunting on stability paraoeterB for flat plate and 10° half -angle cone. Surface- 
to-amoient temperature ratio. 4*0. 



0 2 4 6 6 10 12 14 16 18 20 22 

Fll^t Mach number, Mq 

Figure 7. - CoTt^jarlBon of Inminar-heat-tranaf er ratios for spherical and conical noses. 3urface-to-amhient 
static -temperature ratio, 1.0. 
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